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Abstract 

We study plasma oscillations in a flux tube of the dual superconductor model 
of 't Hooft and Mandelstam. A magnetic condensate is coupled to an electro- 
magnetic field by its dual vector potential, and fixed electric charges set up a 
flux tube. An electrically charged fluid (a quark plasma) flows in the tube and 
screens the fixed charges via plasma oscillations. We investigate both Type I 
and Type II superconductors, with plasma frequencies both above and below 
the threshold for radiation into the Higgs vacuum. We find strong radiation 
of electric flux into the superconductor in all regimes, and argue that this 
invalidates the use of the simplest dual superconductor model for dynamical 
problems. 

I. INTRODUCTION 

The confinement of color in quantum chromodynamics is explained if color electric fields 
form flux tubes. These flux tubes appear readily in the bag model U, but this is a geometric 
picture rather than a dynamical one. 't Hooft [|| and Mandelstam || proposed that, just 
as magnetic flux tubes form in a superconductor 0], a condensation of magnetic charge in 
the QCD vacuum would lead to the formation of electric flux tubes and confinement. This 
is the dual superconductor picture of confinement. 

The analytical study of the confining flux tube depends largely on a classical, Abelian 
model. One can connect this model to QCD with 't Hooft's idea of Abelian dominance ||. 
Starting in QCD, one fixes an Abelian gauge and asserts that the dominant degrees of free- 
dom are Abelian gauge fields (from a Cartan subalgebra) and Abelian magnetic monopoles, 
with weak coupling to the non- Abelian gauge fields which assume the guise of self-coupled 
charged fields. The assumption that the effective interaction of the monopoles causes their 
condensation leads immediately to the dual superconductor picture. The validity of this 
scenario is a subject of research and debate in the lattice gauge theory community ||. 
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In superconductivity, one studies the static structure of magnetic flux tubes via a Landau- 
Ginzburg theory 0. The simplest Landau-Ginzburg Hamiltonian is that of the Abelian 
Higgs model in three dimensions. Classical solutions of the Abelian Higgs theory have been 
considered as models of QCD's electric flux tube as well |]8HT2|. Suganuma, Sasaki, and 
Toki and their collaborators |T3| , |T4|| , using the formalism proposed by Suzuki have 
fixed the effective coupling constants of the Abelian Higgs model by comparing its flux tube 
with phenomenology. There have also been attempts to do so by comparison with the flux 
tube that emerges in lattice gauge theory |I6 |. Success in this program will establish a 
Landau-Ginzburg effective Hamiltonian of QCD. 

In this paper we present a study not of the statics of the confining flux tube, but of its 
dynamics. The flux tube of QCD is not a static object. Creation of flux tubes and their 
subsequent decay through pair creation offer a detailed model of particle creation in e + e _ , 
pp, and pA collisions |17|JT8|| . In the case of nucleus-nucleus collisions, the flux tube has an 
appreciable transverse extent, so that it can be called a "color capacitor" if the color field is 



19| if it is not. qq pairs and gluons are created through 
20|j , and they screen the field in the flux tube while 



coherent across it, or a "color rope" 
the Schwinger pair creation process 
carrying away the energy in the form of hadrons. If the density of created particles is large 
enough, the quarks and gluons form a quark-gluon plasma before the final hadronization 
takes place. 

Field-theoretic analysis of pair creation and back-reaction in the flux tube |2l| has shown 
the buildup of particle density and subsequent plasma oscillations. The methods were ap- 
plied first to a field region of infinite spatial extent, rather than to a finite flux tube. Eisen- 
berg carried out a calculation in a cylindrical flux tube of fixed radius, as suggested by 
the bag model.[] In our study here, the flux tube is a dynamical field configuration, not a 
static geometric object. 

Following |y| and we set up the electric flux tube in classical electrodynamics 
coupled to an Abelian Higgs field via the dual gauge field. Thus the Higgs field represents 
a condensate of magnetic monopoles; the dual Meissner effect confines electric flux to flux 
tubes. Starting from a static flux tube configuration, we release a density of electric charges 
in the system and allow them to accelerate and screen the electric field. The weakening of 
the electric field allows the flux tube to collapse, but the inertia of the charges carries them 
into plasma oscillations that build up the field strength again and force open the tube. The 
coupling of the plasma oscillations to the Higgs field making up the flux tube is the main 
new feature in our work. 

In Section |H] we describe the dual superconductor theory. The model contains an Abelian 
gauge field governed by Maxwell's equations, with coupling to electric charges and currents 
on the one hand and to a magnetic Higgs field on the other. The electric currents are 
treated hydrodynamically; the Higgs field is a classical field that confines electric flux to the 
tube. The coupling to the Higgs field is accomplished by means of a dual vector potential, 
as introduced by Zwanziger p| and used subsequently in || for the static problem and 



This calculation imposed superconducting boundary conditions at the tube surface, instead of 
the more appropriate dual superconductor. 



2 



in |T3| , |14|1 for phenomenological modeling. We impose cylindrical symmetry and we assume 
^-independence along the flux tube in the region far from its ends. Nevertheless, we are 
forced to pay some attention to what happens near the ends of the tube in order to define 
potentials unambiguously. 

We present numerical results in Section [[II]. We set the parameters of the theory in 
both the Type I and Type II regimes, and show how plasma oscillations arise. We consider 
plasma frequencies both above and below the cutoff for radiation into the Higgs vacuum. 
The plasma oscillations are accompanied by changes in the radius of the flux tube as the 
electric pressure from within decreases through screening, only to increase again as the 
currents overshoot. 

We believe that this paper is the first to address the dynamics of an electric flux tube 
within this model.0 It is amusing to note that this physical situation has no counterpart in 
superconductivity. While our magnetic Higgs field appears in the usual Landau- Ginzburg 
theory as an electrically charged condensate, our electric currents cannot appear there be- 
cause there are no magnetic monopoles in nature. An Abrikosov flux tube has to run all 
the way to the boundary of the sample, where it joins onto the external magnetic field. Our 
electric flux tube, on the other hand, has a finite (but large) length, and the charges at its 
ends can be screened by the electric currents that flow in it. We review in the appendix 
the well-known phenomenon of flux quantization and why it has no effect on a flux tube of 
finite length. 



II. FIELD EQUATIONS IN CYLINDRICAL GEOMETRY 

Our classical model for the flux tube is electrodynamics coupled to a scalar magnetic 
monopole field. The equations of motion are Maxwell's equations with both electric and 
magnetic sources, plus the Klein-Gordon equation for the monopoles. The latter is coupled 
not to the usual vector potential A M but to the dual potential J23| B^, and contains a 
self-interaction that puts it into the Higgs phase. 

We make a cylindrically symmetric ansatz for the fields, and furthermore assume z- 
independence far from the sources at the ends of the tube. We put the electrically charged 
sources on the z axis at z = ±z , an d we assume that all the space-charge effects are localized 
there. The central region is defined by \z\ < z\, where z\ < z is chosen to exclude the end 
regions (see Fig. [I]). 



A. Maxwell's equations 

Maxwell's equations with electric and magnetic currents are 

d,F^=f e (2.1) 
d,F^ = f g , (2.2) 



Loh et al. |24|] have simulated breaking of the flux tube in the Friedberg-Lee model [25 which 
posits a confinement mechanism unrelated to monopole condensation. 



where F^ v = ±^ uXa F Xcr , Eq. ( |2.2|) replaces the Bianchi identity of ordinary electrodynamics. 
We can solve it with a vector potential, but there is a new term, 

pi* = qh A » _ d v A » + e ^G Xa . (2.3) 

Defining an arbitrary vector n M , one finds that ( |2.2| ) is solved if 

GT = -n tx {n-d)- x f g , (2.4) 

where (n-<9) _1 represents an integration from infinity with suitable boundary conditions. The 
field equation for then follows from Q2.1|) . Solving Maxwell's equations in the opposite 
order, the same field strength can be represented alternatively by a dual vector potential, 
according to 

piw = d n B v _ ff> B n + e ^Mxc , (2.5) 

with 

M^ v — — n M (n - d)~ x jl , (2.6) 



in order to satisfy ( |2.1| ). Now the field equation for follows from ( |2.2| ). We will return to 
the vector potentials below. 

In three-dimensional notation, 

V • E = p e (2.7a) 

dE 

VxH- — =j e (2.7b) 

V ■ H = Pg (2.7c) 



V x E + — = -j s (2.7d) 



In the aforementioned central region \z\ < z%, we make the ansatz 

E = E{r,t)z (2.8a) 

H = H(r,t)9 (2.8b) 

j e = j e (r,t)z (2.8c) 

h=j 9 (r,t)§ (2.8d) 
Ampere's Law ( |2.7b| ) has only a z component, 

4 



and Faraday's Law ( |2.7d|) has only a component, 

BE OH , 

The two Gauss laws show that the charge densities are zero, 

Pe=Pg=0. (2.11) 

Eqs. (|2.9| ) and ( 2.10|) are the equations of motion for and if. 

Near the ends of the flux tube, \z\ > z±, we must allow for new components E r and j er , 
as well as for p e , and all fields will be ^-dependent. 



B. Vector potentials 

In order to calculate the time evolution of the matter fields, we need the vector potentials. 
In fact, we need only the dual potential fl26 |, which enters the field equation for the 
monopole field. From (|2.5|) and 



E 
H 



Iz-VHzxje) 



(2.12) 
(2.13) 



-V x B + z(z- V)~V 

— f 

where we have chosen = (0, z). We choose the gauge Bq = 0. Inverting ( |2.12| ) gives 

B = (z-Vr 1 (zxE). (2.14) 
Our ansatz for E then gives B = B 9, with 



B(r,z) = f E r (r,z')dz' 

J — oo 



(2.15) 



This integral only gets contributions from the regions around the sources, since E r = far 
from the charges. This makes B independent of z in the central region \z\ < Z\. We can 
relate the integral in ( |2.15| ) to the charge distribution near the sources as follows. We take a 
cylindrical surface of radius r with ends at — oo and —z\. Gauss' Law gives for this surface 



-21 



Q(r) = 1-iir \ E r (r, z) dz + 2vr / E z (r', -Zi) r' dr' 



(2.16) 



where Q(r) is the charge inside the cylinder (composed of the original source plus the space 
charge around it). Noting that E z (r', —Zi) is the ^-independent E z (r'), we obtain the relation 



and its inverse 



Br) 



E = E, 



2-rrr 



- f EJr')r'dr' 
r Jo 

15 1 dQ 

-- — {rBj H — . 

r or 2nr or 



(2.18) 



Given E and some model for the charge Q(r) (see below), we can calculate B. 
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C. Magnetic monopoles 

The magnetic monopoles are represented by a classical Klein-Gordon equation with a 
Higgs potential, 



where 



D^D^^ + X(\^\ 2 -v 2 )^ = 0, 



D^ = d^-igB,. 



We make the ansatz ip = p(r,t)e lX9 ^ 9 \ giving 



OH- _ -— r — - (-— - igB 

dt 2 r dr dr \r dO 



pe ix + X (p 2 - v 2 ) pe 



(2.19) 



(2.20) 



(2.21) 



We assume Xg — n @ (f° r a A ux tube with n units of electric flux), so 



d 2 



1 d d 



dt 2 r dr dr 



-gB^j 2 p+ \(p 2 -v 2 )p = 



(2.22) 



This is the equation of motion for p. It feeds back into Maxwell's equations via the magnetic 
current, 



• h = 2gp 2 (V Xg ~ gB) 
which contains only a 9 component, 

j g = 2gp 2 (- - gB 



(2.23) 



(2.24) 



As shown in the Appendix, n represents electric flux coming in from infinity, so that we 
will set n = for simplicity. The flux generated by Q at the ends of the flux tube is not 
quantized. 



D. Charged matter 

We represent the electrically charged matter by classical two-fluid magnetohydrodynam- 
ics. The positively and negatively charged fluids both have particle density n e (r,t) (so that 
p e = 0) and their velocities are v ± . The fluid motion is determined by Euler's equations 



m 



~dt 



+ (v± • V)v J 



±eE ± ev ± x H — —VP 

n e 



(2.25) 



and the continuity equation 



dt 



+ V • (n e v ± ) = 



(2.26) 
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The electric current is given by 

(v + - v-) . (2.27) 



We assume in the central region that all quantities depend only on r and t, so that 
±d_ 

r dr' 



v 1 * 1 • V = vf-K-. Recalling that E = E z and H = H 9, we have 



and 



9 ± ± d ± ±u IdP 

m—vt = —mvt—vt =f evtn — 

dt or n e or 

d d 

m—vf = -mvf—vf ±eE± evfH (2.28) 
Ot Or 



°"' ~vt^ ~ ^Urvf) ■ (2.29) 



dt r dr r dr 

Not surprisingly, these equations allow the ansatz 



Vt = V„ = V T 



we find 



d d TT 1 dP 

\—v r = -mv r —v r - ev z H — 

ot or n P Or 



9 d 
m—v z = —mv r —v z + eE + ev r H 
ot or 

d Jl± = - Vr ^-!hl( rVr) . (2.30) 
dt r dr r dr r 

These equations feed back into Maxwell's equations via the current, which is in the z direc- 
tion and has the magnitude 

j e = 2n e ev z . (2.31) 

We need an equation of state to relate P to n e . Having in mind a quark-gluon plasma, 
we choose the equation of state of a relativistic ideal gas with the appropriate number of 
massless fermions and bosons. We set the chemical potential to zero and write all quantities 
in terms of the temperature, 

^ 2 4 

P = 37 — T 



90 

n e = 37^T 3 . (2.32) 



7T 

Then eqs. (|2.30| ) can be written in terms of the velocities and temperature as 
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d d 4vr 4 dT 

m—v r = —mv r —v r — ev z H rr^^r- 

dt dr 90C(3) dr 

9 d 
m—v z = —mv r —v z + eE + ev r H 
dt or 

dT dT T d 

dt r dr 3r dr r 

The advantage of a hydro dynamic description of the charged matter is that the fluid 
interacts directly with the field strengths E and H. If we were to consider classical or 
quantum field theory instead of hydrodynamics, we would need the vector potential A. In 
our geometry, the simplest ansatz would require two components, A r and A z . The potential 
would necessarily be z-dependent in the central region; moreover, the simplified treatment 
of the ends of the flux tube (see below) would be impossible. 



E. The ends of the flux tube 

For study of the central region, all we need to know about the complex regions at the 
ends of the flux tube is Q(r), which represents the charge density built up near — zq by the 
current j e . (The region near z = +zq is, of course, a mirror image of the region near —Zq.) 
Without dealing in detail with the motion of charges near the ends of the flux tube, we can 
guess at a few models: 

1. Point charge: Here we just assume that all the current merely accumulates in a point 
charge at z — —zq. Then Q(r,t) is independent of r, and charge conservation gives 

dO f°° 

= -2tt / j e (r)rdr. (2.34) 
dt Jo 

2. Surface charge density: Here we let the charge pile up on a plate near z = —z , giving 
a surface charge density cr(r,t). Then 



Q(r,t) = 2ti I a(r',t)r'dr' (2.35) 

and 



o 



!=-*■ (2 - 36 » 

Initial conditions have to be specified for o to represent the original source of the flux 
tube. 

3. Both: Keeping the initial charge pointlike, we put only the accumulated charge into 
a. Then 

Q(r, t)=Q + 2vr f a(r', t) r' dr' (2.37) 
Jo 

and a is determined by ( |2.36| ) as above. The initial condition is a = 0. Note that in 
this case the space charge will never exactly cancel Qq and thus any plasma oscillations 
will be asymmetric. 
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The initial conditions must in any case satisfy Gauss' Law, 

roo 

2tt / E(r)rdr = Q(oo) , 
Jo 

if there is no flux coming from infinity (i.e., n = 0). 



(2.38) 



F. Initial conditions 



We start off the system in the configuration of a static flux tube with a stationary fluid 
in it. We combine the static limit of Q2.10 ) with ( 2.18Q and ( p. 24 ) to give 



d 1 d 
dr r dr 



rB(r) - —Q(r) 

Z7T 



2g 2 p 2 B 



[2.39) 



where we have chosen n = 0. Eq. ( p. 171 ) gives the boundary condition B(0) =0 unless Q 
includes a point charge Qq; in that case 



B(r) 



2nr 



(2.40) 



as r — > 0. Eq. ( |2.17p together with ( |2.38| ) gives rB(r) — > at infinity. 



The static limit of the Klein-Gordon equation ( 2.22 ) is 

1 d d 



--^r— + g 2 B 2 
r or or 



p + A (p 2 - v 2 ) p = 



(2.41) 



with the boundary conditions that p is zero at the origin and tends to p = v at infinity. 
Focusing on the case where Q is a point charge, we define the reduced field b(r) via 



giving 



2nr 



b" 2g 2 p 2 b = 



(2.42) 



(2.43) 



with 6(0) = 1 and 6(oo) = 0, and 



- P 



2n r 



+ X(p 2 -v 2 



p = 



[2.U) 



These equations determine p and B, and hence E, for the static flux tube. Clearly H = 
here. 
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III. PLASMA OSCILLATIONS 



We determine the time evolution of the system through the following system of equations. 
The Maxwell equations ( |2.9| ) and fl2.10| ) give E and H; the MHD equations (|2.33|) give v and 
T, and hence j e via fl2.31|) and ( p.32|) . Eq. ( |2.34|) gives the charge Q whence ( |2.17|) gives the 
vector potential B. The scalar field p evolves according to the Klein-Gordon equation Q2.22|) , 
and (p. 24 ) gives the magnetic current j 9 , the last ingredient for the Maxwell equations. The 
initial conditions, as noted, consist of the static flux tube with an initial value of Q, and of 
a static fluid distribution specified by T{r). 

We choose four sets of parameters, listed in Table |. The magnetic charge g and the 
vacuum expectation value v = p(r — > oo) determine the vector mass my = \[2gv and its 
reciprocal, the London penetration depth A^ = m v . The scalar self-coupling A determines, 
along with v, the (dual) Higgs mass m^ = \/2Xv and the Ginzburg-Landau coherence length 
£ = m^ 1 . The ratio of these lengths is k = = vA/S'- If this value is smaller than one 
then the superconductor is of Type I; otherwise it is of Type II. 

In ordinary superconductors, flux tubes are observed only in Type II materials, when the 
applied magnetic field lies between the critical values H c \ and H C 2 and penetrates via creation 
of an Abrikosov lattice. Type I materials, on the other hand, expel the field entirely as long 
as superconductivity persists. The situation would be different if one were to introduce 
a pair of magnetic monopoles. Then the magnetic field would be forced into a flux tube 
between the monopoles, even in a Type I material. The difference between Type I and II 
would lie in the stability of the flux tube against splitting into a lattice (a "rope" ) of smaller 
tubes: Type I tubes would be stable, while Type II tubes would split. This splitting is 
presumably inaccessible from our cylindrically symmetric ansatz. 

The values of A and g fixed phenomenologically in |TJ|] are in the Type II region, and we 
choose these values for our Type II cases. For Type I, we increase g and decrease A to make 
k < 1. We take v = 126 MeV from []13| as well. We fix the electric charge e according to the 
Dirac quantization condition eg/An = 1, even though this is not meaningful for continuum 
hydrodynamics. 

The charged fluid presents a conundrum. Our intuition about confinement suggests that 
the superconducting vacuum should expel this matter, and thus confine it to the flux tube. 
(This happens in the Friedberg-Lee Model 125].) In this Abelian theory, however, the locally 
neutral fluid is not confined. If the initial conditions contain a fluid inside the flux tube 
only, it will flow outward in a hydrodynamic rarefaction wave. Perhaps the fluid outside the 
flux tube may be interpreted as a hadronic fluid. In any case, we prefer not to superimpose 
radial hydrodynamic flow on the oscillations of the flux tube. Thus we choose a uniform 
initial fluid density, outside the flux tube as well as inside. Having in mind a quark-gluon 
plasma, we set the initial temperature to be spatially uniform with To = 150 MeV. 

The initial value Qq of the charge at the ends of the flux tube gives the initial amount 
of electric flux. In a pp collision, one would have Qq ~ e, corresponding to the creation of 
a flux tube in the fundamental representation; an AA collision between heavy nuclei would 
give |19|>23 Qo ~ A 1 ' 3 e, which reaches Q ~ 6e for large nuclei. A large value of Q will give 



a thick flux tube and large-amplitude plasma oscillations, and thus enhance the non-linear 
effects due to the Higgs coupling. Unfortunately, very large values of Qo are outside the 
range of stability of our numerics; we choose values that are as large as possible given this 
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constraint. 

The plasma frequency is 




UJ v = \—— , 



where the factor of two reflects the fact that we have two fluids. For given n e and e, we fix 
m so as to tune u p to either side of the vector mass my; thus the plasma oscillations will 
occur at frequencies either above or below the threshold for radiation. We expect a priori 
that m will be in the neighborhood of the dynamical mass eT of light particles in a heat 
bath. The parameters used in [|13[] give a rather small vector mass, and with e = 47r/g the 



condition u p < my leads to an unreasonably large value of m. Instead we choose in this 
case to lower the value of the electric charge to e = 1. 

Case 1. Fig. ||] shows plasma oscillations in the on-axis electric field for the Type I 
superconductor with u p < my. The oscillations are clearly nonlinear, with varying amplitude 
and misshapen waveform. Fig. || contains snapshots taken at the times of maxima in E(r = 
0). These snapshots show radiation emanating from the flux tube in both E and p. The plots 
of rE show that the oscillating electric flux in the outgoing wave is as large in amplitude 
as the initial flux distribution. It is difficult to tell, however, whether the wave amplitude 
decays as r -1 / 2 , as expected for a true propagating wave, or as r _1 , which would make the 
wave evanescent. 

Since the frequency of the plasma oscillations is below the vector mass, this cannot 
be simple linear radiation. The third row of fig. |3], indeed, shows strong coupling of the 
nonlinear p field to the electromagnetic wave. The origin of this coupling, shown in fig. [|, is 
the collapse of the flux tube under the pressure of the p field when the electric field is weak. 
The figure shows snapshots of p taken in the first half-cycle of the evolution. The middle 
snapshot, taken when E(r = 0) = 0, shows appreciable narrowing of the flux tube core even 
as p is driven downward at larger radii in the first oscillation of the outgoing wave. 

Case 2. Oscillations in the Type I superconductor with a larger plasma frequency, 
u p > my, are shown in figs. f| and || Here the irregularity of the waveform of E(r = 0) 
is due to ordinary radiation in interaction with a more-or-less static flux tube wall. Again 
there is appreciable electric flux radiated outward, but there is little effect on the p field. 
This is because the latter cannot respond when driven by the high-frequency electromagnetic 
oscillations. We note that the larger value of gQo here creates a thicker flux tube and keeps 
more flux in the flux tube despite the radiation. 

Case 3. The Type II superconductor with u p < my is similar in behavior to the Type 
I system. The smaller plasma frequency makes it difficult to calculate through more than 
two or three oscillations, in spite of the high stability of the Crank-Nicholson algorithm we 
use. In fig. [7| we see strong deformation of the waveform. Fig. |8| shows radiation similar to 
the Type I case, but the perturbation of the p field is less apparent. 

Case 4- Finally, the corresponding plots for the Type II superconductor with u v > my 
show considerable structure. Recall that this parameter set is that used in |13| for phe- 
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nomenological fits to static quantities.0 The oscillations shown in fig. |9] are too rapid to 
couple strongly to p, and we see strong electromagnetic radiation (figs. 10 and 11) accom- 
panied by weak high-momentum disturbances in p. 



IV. DISCUSSION 



Our numerical results raise difficulties for the dual superconductor picture of confinement. 
As a static model, the dual superconductor does indeed form flux tubes that confine charges 
in string-like configurations. Once the dynamics are examined, however, the lack of absolute 
color confinement becomes apparent. While the motion of neutral particle matter outside 
the flux tube may be passed off as the emission of color-neutral hadrons, the radiation of 
appreciable electric flux cannot. The electric field in 't Hooft's Abelian projection is after 
all a color field, representing a coherent, colored gluon state. 

It was predictable that oscillations with u p > my would radiate into the Higgs vacuum, 
since the photon does have a less-than-infinite mass. One might be tempted to restrict 
application of the dual superconductor model to situations where frequencies are much less 
than my. We have seen, however, that this is insufficient. Even low- frequency plasma 
oscillations, where radiation should be impossible, succeed in spreading out the electric flux 
via nonlinear effects. Thus the model is not reliable even in this regime. It is possible 
that the model can be made applicable to low-frequency physics by choosing a monopole 
Lagrangian more general than the ip A of the simplest Landau- Ginzburg theory. Absolute 
confinement, however, will never be realized in this way. (The introduction of a strongly 
nonlinear dielectric constant is what enables the Friedberg-Lee model to confine all color 
fields, but this model has not been derived from QCD.) 

't Hooft's Abelian reduction rests on the identification of the important degrees of free- 
dom in an Abelian gauge. It is supposed that the magnetic monopoles have a strong 
self-interaction, leading to their condensation; that the Abelian gluons, belonging to the 
Cartan subalgebra, turn this condensate into a superconductor; and, most important, that 
the off-diagonal gluons are irrelevant to confinement, except insofar as they screen zero- 
triality states. These off-diagonal gluons, however, retain all the self-couplings of the orig- 
inal non-Abelian gauge theory. We conjecture that they are essential to understanding 
time-dependent phenomena related to confinement .[] 
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[ PI , and hence no plasma frequency. 
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APPENDIX: FLUX QUANTIZATION 

Flux quantization, as it turns out, affects only flux coming in from infinity and not the 
flux due to the charge Q at the ends of the flux tube. The quantization condition comes from 
demanding that the total energy be finite ||28|. The energy of the monopole field contains 



the term 

2 



^static = I d X 



T) B i) . (Al) 



As x 2 + y 2 — > oo, we have ip ~^ ipoe lx ^ and B — > B 9, and thus the integrand approaches 

eoo^o 2 -(^-<^) 2 • (A2) 
In order that J°° rdrdOe be finite, we must have at large r that 

§ = (A3) 

Since \ is only determined up to a multiple of 2tc, it can gain such a multiple when one goes 
around the circle. Integrating (|A3|) over the circle at fixed r gives 



2vm = g J BrdO = g j>& ■ dl , (A4) 

where n is an integer. In view of ( |2.12| ), 

j> B • d\ = y (V x B) • dS = - J E • dS + Q = -$ B + Q (A5) 

where $^ is the total electric flux. Thus 

§ E = Q - 27in/g . (A6) 

Only the external flux, the flux that does not end at Q, is quantized in units of 27r/g. 
Choosing = Q, i.e., no external flux, means setting n = 0. 
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TABLES 



TABLE I. Parameters used in numerical simulations, v = 126 MeV in all cases, and the initial 
temperature is Tq = 150 MeV. 





Tvue I 


Tvoe I 


Type II 


Type II 






LO^ ^> TTi s \r 

- — ' ' c y 




LJ^, "> 777,1/ 


A 


10 


10 


25 


25 


a 

it 


5.5 


5.5 


2.3 


2.3 


e 


2.3 


2.3 


1.0 


5.5 


my (MeV) 


975 


975 


411 


411 


m H (MeV) 


565 


565 


893 


893 


Al (fm) 


0.20 


0.20 


0.48 


0.48 


£ (fm) 


0.35 


0.35 


0.22 


0.22 


K 


0.58 


0.58 


2.18 


2.18 


m (MeV) 


250 


50 


250 


50 


uj p (MeV) 


790 


1766 


343 


4192 


T p = 2tt/lo p (fm) 


1.57 


0.70 


3.60 


0.30 


Qo 


2.3 


4.6 


1.0 


4.6 
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FIGURES 
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FIG. 1. Geometry of the flux tube. 
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FIG. 2. Electric field on the z axis, for Type I superconductor with uj p < my. 
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Type I, to < m v 



t = 0.0 t = 3.0 



t = 7.8 




FIG. 3. Snapshots of the electric field E and monopole field p at times of maxima in the on-axis 
field E(r = 0): Type I superconductor, lo p < my- 
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p(r,t):Type I, co < m v 
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FIG. 4. Snapshots of p(r, t) in the first oscillation: Type I superconductor, uj p < my- 



2.0 



1.0 



0.0 



-1.0 



-2.0 



0.0 



E(r=0,t) 

Type I, co p > m v 



2.0 4.0 6.0 

t (fm/c) 



8.0 



FIG. 5. Electric field on the z axis, for Type I superconductor with lj p > my. 
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FIG. 6. Snapshots of the electric field E and monopole field p at times of maxima in the on-axis 
field E(r = 0): Type I superconductor, lo p > my. 
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E(r=0,t) 



Type II, co p < m v 




0.0 2.0 4.0 6.0 8.0 
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FIG. 7. Electric field on the z axis, for Type II superconductor with uj v < my- 
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Type II, co p < m v 

t = 0.0 t = 3.4 t = 7.0 




FIG. 8. Snapshots of the electric field E and monopole field p at times of maxima in the on-axis 
field E(r = 0): Type II superconductor, lj p < my. 
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FIG. 9. Electric field on the z axis, for Type II superconductor with uj p > my- 
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Type II, co p > m v 

t = 0.0 t= 1.1 t = 2.2 




FIG. 10. Snapshots of the electric field E and monopole field p at times of maxima in the 
on-axis field E(r = 0): Type II superconductor, oj p > my. 
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Type II, co p > m v 
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